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ABSTRACT 

We have analysed 5-epoch GOODS HST-ACS B, V, I775 and z datasets (VI. release), 
in conjunction with existing VLT-ISAAC imaging in the J, H and Kg bands, to derive 
estimated redshifts for the sample of 198 Extremely Red Galaxies (ERGs) with Kg < 
22 (Vega) and (I775 - K^) > 3.92 selected by Roche et al. (2003) from 50.4 arcmin^ 
of the GOODS/CDFS field. We find that, at this depth, the ERG population spans 
\ the redshift range 0.5 < Zphot < 4.75 and over two decades in mass (~ 3 x IO^Mq 

^ ■ to ~ 3 X lO^^M©). Our results show that the dust-corrected red envelope of galaxy 

\ evolution is well modelled by a starburst at redshift Zf = 5 followed thereafter by 

tJ" . passive evolution. We explore the evolution of the ERG luminosity function (LF) from 

' redshifts (zphot) = 1-0 to (zphot) = 2.5 and compare it with the global Kg-band LF at 

redshifts 1 < Zphot < 2. We find that the bright end of the ERG LF does not decrease 
from redshifts (zphot) — 2.0 to (zphot) = 2.5 and we connect this fact with the presence 
of progenitors of the local L > L* population at redshifts Zphot > 2. We determine 
lower limits of = (6.1 ± 1.9) x lO^^Mpc"^ and pc ^ (2.1 ± 1.1) x 10^^ Mpc'^ 
1-^ ' on the comoving densities of progenitors of local massive galaxies already assembled 

at redshifts (zphot) = 2.5 and (zphot) = 3.5, respectively. We have investigated the 
Q ■ existence of high-redshift Lyman break galaxies massive enough to be included in 

$H ' this ERG sample. Out of an initial list of 12 potential very high redshift candidates, 

lyj , we have identified 2 ERGs which have a high probability of lying at Zphot > 4. We 

discuss the advantages of multi-colour to single-colour selection techniques in obtaining 
reliable lists of very high-redshift candidate sources, and present revised lower redshift 
estimates for sources previously claimed as potential z > 5 dropouts in recent studies. 

■ Key words: galaxies: evolution - galaxies: formation - galaxies: high-redshift 



1 INTRODUCTION 

The study of 'extremely red galaxies' (ERGs) is of im- 
portance for setting constraints on the first epoch of mas- 
sive galaxy formation, for exploring the population of dust- 
reddened starbursts, and for testing models of galaxy forma- 
tion in general. Usually defined as objects with R — K > 5 
or I — K > 4 (Vega), these galaxies have been widely studied 
since their discovery by Elston, Rieke & Rieke (1988). The 
traditional picture identifies the ERGs either with passively 
evolving elliptical galaxies, or with young starbursts strongly 
reddened by dust. This mixture has been confirmed by spec- 
troscopic (Dunlop et al. 1996; Soifer et al.l999; Cimatti et 
al. 2002; Saracco et al. 2003) and sub-millimetre surveys 
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(Cimatti et al. 1998; Dey et al. 1999; Wehner, Barger & 
Kneib 2002). However, the latest photometric and morpho- 
logical studies of these objects indicate that the ERG popu- 
lation is perhaps more complex than previously considered 
(Small et al. 2002; Miyazaki et al. 2002; Van & Thompson 
2003; Moustakas et al. 2004). 

Most previous samples of ERGs have been selected to a 
limiting magnitude of Kg = 20 or brighter, and their inferred 
redshift distributions peak at z ~ 1 — 2 . New deeper obser- 
vations reaching Kg = 22 and beyond (Maihara et al. 2001, 
Saracco et al. 2001) are uncovering ERGs at higher redshifts. 
In particular, a population of galaxies with very red near- 
infrared colours J — Ks > 3 has been discovered to lie at red- 
shifts z > 2 (Totani et al. 2001). Saracco et al.(2004) claim 
that three massive galaxies of this kind in the Hubble Deep 
Field South (HDFS) at (zphot) = 2.7 account for about 40% 
of the comoving density of the local early- type L > L* galax- 
ies. However, these studies are based on small area surveys 
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Figure 1. Schematic diagram of tlie sub-region covered by 
deep VLT/ISAAC near-infrared imaging (sliaded area) within the 
larger GOODS/CDFS field with full HST-ACS coverage. 



also demonstrated that the strong clustering displayed by 
the ERG population extends to the faintest Ks magnitudes. 
In this work, we present photometric redshift estimates for 
the ERGs in the Roche et al. sample and explore the result- 
ing implications for the nature and evolution of ERGs. 

The layout of this paper is as follows. First, in Sec- 
tion 2, we briefly review the sample selection and explain 
the details of the multiwavelength photometry. In Section 
3 we discuss the photometric redshift techniques applied, 
and present examples of individual redshift determinations. 
Then, in Section 4 we present our results and discuss raw 
and dust-corrected Hubble diagrams (Kg — Zphot), the red 
envelope of galaxy evolution, the evolution of the ERG lu- 
minosity function (LF), and derived lower limits on the co- 
moving densities of progenitors of local L > L* galaxies. We 
discuss the individual properties of each of our high-redshift 
(zphot > 4) candidate sources and revise the nature of those 
included in the Hubble Ultra Deep Field (UDF), in the light 
of the new publicly available ultra-deep images taken with 
the HST-ACS and the HST-Near Infrared Camera and Multi 
Object Spectrometer (NICMOS). We also review the exis- 
tence of an ERG cluster reported by Roche et al. (2003). Fi- 
nally, we present some concluding remarks in Section 5. We 
adopt throughout a cosmology with Ho = 70 kms~^Mpc~^, 
r^M = 0.3 and = 0.7. 



which cover a few square arcmin of the sky. Deep Ks-band 
observations in wider areas are required to obtain significant 
samples of high redshift ERGs and to determine more accu- 
rate lower limits for the fraction of massive galaxies already 
assembled at different redshifts. 

The Great Observatories Origins Deep Survey 
(GOODS) (Dickinson et al. 2003) is providing unprece- 
dented multiwavelength data in ~ 320 arcmin'^ centred 
on the Chandra Deep Field South (CDFS) and Hubble 
Deep Field North (HDFN). Within the GOODS/CDFS 
field, Roche, Dunlop & Almaini (2003) selected 198 ERGs 
with Ks < 22 (Vega) and I775 - > 3.92 from the 50.4 
arcmin^ for which deep near infra-red (near-IR) data have 
been obtained with the Infrared Spectrometer and Array 
Camera (ISAAC) on the 'Antu' Very Large Telescope 
(Antu-VLT). This is the deepest significant sample of ERGs 
selected to date, and is the subject of the present study. 
The region covered by the ISAAC observations within 
the GOODS/CDFS field is shown in figUl Other previous 
studies of ERGs reaching Ks ~ 20 exist in the CDFS. 
For example, the K20 survey has yielded spectroscopic 
redshifts for ERGs in 32.2 arcmin^ of this field (Cimatti 
et al. 2002). More recently, the GOODS team (Moustakas 
et al. 2004) selected a sample of ERGs from 163 arcmin^ 
of the GOODS/CDFS field with data from the Advanced 
Camera for Surveys (ACS) on board the Hubble Space 
Telescope (HST). The area covered by this brighter ERG 
sample is three times greater than that covered by the 
Roche et al. sample studied here. However, the Moustakas 
et al. sample only reaches a depth of Kg ~ 20.2 (equivalent 
to Ks — 22, AB), and so the present study probes a different 
(and complementary) region of parameter space to that 
already explored by the GOODS team. 

Roche et al. (2003) studied the number density of ERGs 
within their Ks = 22 sample and compared the results with 
the predictions of different galaxy formation models. They 



2 THE SAMPLE-MULTIWAVELENGTH 
PHOTOMETRY 

GOODS observations include optical and near-IR imaging 
in the B, V, I775 and z bands with the ACS/HST and 
J, H and Ks bands with the ISAAC- VLT. Roche et al. 
(2003) used ISAAC K^-band data and ACS 1775-band data 
to select a sample of 198 ERGs with Ks < 22 (Vega) and 
I775 — Ks > 3.92. They performed the source extraction in 
the Ks band and photometry in the I775, J, H and Ks bands 
using the public code 'SExtractor' (Berlin & Arnouts, 1996). 
All the selected ERGs are at least 3a detections on the Ks- 
band images. Further details of their photometric measure- 
ments from the ground-based data and their preliminary 
measurements on the 1775-band ACS images are given in 
their paper. The aperture magnitudes measured on ISAAC 
images by Roche et al. (2003) and used here correspond to 
2"-diameter in all cases. 

We performed photometric measurements for the four 
ACS bands on the stacked GOODS 5-epoch images (VI. 
HST/ ACS data release), using the corresponding weighting 
maps. We looked for counterparts of the Ks-selected ERGs 
on the I775 and z-band images within an angular radius of 
1". For B and V bands we restricted the search of counter- 
parts to 0.5" to minimise the presence of interlopers. We 
also performed the magnitude measurements using the pub- 
lic code 'SExtractor'. Many ERGs appeared as non-detected 
in the B and V bands, and some of them as non-detected 
in the I775 or z bands. The algorithms for computing pho- 
tometric redshifts depend critically on the adopted value 
of limiting magnitudes in the case of non detections, and a 
careful treatment of potential 'dropouts' is necessary to pre- 
vent the photometric redshift algorithms from finding fake 
high redshift sources. Therefore, in each case of a 'SExtrac- 
tor' non-detection, we used the IRAF task 'phot' to mea- 
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Figure 2. Example spectral energy distributions (left-hand plots) and redshift probability density distributions (right-hand plots) for 
three ERGs in the Roche et al. sample, a) an ERG with 'hyperz' primary solution Zphot = 1-40, b) an ERG with primary solution 
2iphot = 2.52, c) an ERG with primary solution Zphot = 3.10. The circles in the SED plots correspond to the measured aperture 
magnitudes in each filter. The solid-line curves indicate the best-fit template in each case. 
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sure aperture magnitudes centred at the ERG Ks-band po- 
sition. We measured 2"-diameter aperture magnitudes on 
each ACS filter for those ERGs with zj" < 26 (AB). Beyond 
Z2" — 26 (AB), we found that little flux was missed using 
1" instead of 2"-diameter aperture magnitudes. Thus, for 
objects with Z2" > 26 (AB), we used l"-diameter aperture 
magnitudes to reduce the random errors, and corrected the 
values systematically using the average offset between the 
2" and l"-diameter aperture magnitudes measured in each 
A = 1 magnitude bin. After this procedure, only a few ob- 
jects remained as formally non-detected in the V or redder 
ACS bands. 

The photometry of three ERGs included in an initial 
list of potential very high redshift sources has been revised 
on the new available ACS and NICMOS ultra deep optical 
and near-IR images of the Hubble UDF (cf. Section 4.6). 
The procedure we followed to perform photometric mea- 
surements on these images is similar to the one used for 
the 5-epoch stacked GOODS HST/ACS images. 



3 REDSHIFT ESTIMATIONS 

We computed photometric redshifts for the ERGs in the 
Roche et al. (2003) sample with the public code 'hyperz' 
(Bolzonella, Miralles & Pello, 2000), using the seven pass- 
bands described in Section 2 (B, V, I775, z, J, H and Ks) 
and the GISSEL98 hbrary of Bruzual & Chariot (1993). The 
code 'hyperz' computes the probability of an object being 
at a given redshift Zphot by searching for the best (x^) fit 
to its photometric spectral energy distribution (SED) pro- 
vided by any of the template spectra available in the library. 
The redshift which yields the highest probability is known 
as the 'hyperz' primary solution and, in the following, we 
refer to the best-fit SED as the one corresponding to this 
solution. Although this primary solution is the one used in 
most cases, 'hyperz' also allows one to construct a proba- 
bility density distribution in redshift space for each object. 
This is very useful for revealing possible degeneracies in the 
redshift determination, and for assessing the significance of 
the primary solution. To account for dust obscuration, we 
applied a Calzetti et al. (2000) reddening law, allowing the 
extinction in the V-band (Av) to vary between and 1. For 
the reddest objects, we ran 'hyperz' again allowing Av to 
vary between and 3. This did not produce any substan- 
tial changes in the redshift estimations, but does have some 
impact on the estimated absolute magnitudes and ages. 

We also used the public 'Bayesian photometric red- 
shifts' (BPZ) code produced by Bemtez (2000) to obtain a 
second, independent set of redshift estimates for the ERGs 
in our sample. This code computes photometric redshifts 
using a Bayesian approach (and also has the option of us- 
ing a maximum likelihood technique like the one used by 
'hyperz'). The agreement between the redshifts estimated 
by 'hyperz' and the BPZ code is very good in most cases. 
Accordingly, with a few exceptions discussed individually in 
later sections, we have adopted the 'hyperz' results for the 
ERG sample simply because it provides not only the redshift 
probability distribution, but also several other useful quan- 
tities in the output (e.g. k-corrected absolute magnitudes in 
a selected filter, V-band dust extinction and age). 

Figure|5|shows examples of the best-fit SEDs (left-hand 
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Figure 3. Redshift distribution of the ERGs in Roche et al. 
sample: a) taking into account only 'hyperz' primary solutions, 
b) taking into account the probability density distribution of each 
ERG in the redshift space and small corrections due to the in- 
completeness of the sample. In the upper plot a), the redshift 
distribution of the shallower sample of ERGs selected by Mous- 
takas et al. (2004) in a wider area of the GOODS/CDFS has been 
added for comparison (shaded histogram) , after division by a fac- 
tor of 3.08 in order to normalise their distribution to produce the 
same number of objects in the redshift bin centred at Zpjj^t = 1 
as found in the present study. 



panels) and redshift probability density distributions (right- 
hand panels) for three of the ERGs in our sample. In the 
SED plots, the circles correspond to the measured aperture 
magnitudes in each filter. The solid-line curves indicate the 
best-fit templates. The SED shown in fig. |2l-a) corresponds 
to an object with 'hyperz' primary solution Zphot ~ 1.40. 
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The accuracy in the fitting and the small error bars in the 
photometry are reflected in a relatively small uncertainty in 
the redshift estimation. The SED shown in fig. |3b) corre- 
sponds to an object with primary solution Zphot = 2.52. In 
this case, the larger error bars in the photometry produce a 
wider redshift probability density distribution. Finally, the 
SED in fig. |5J-c) corresponds to an ERG with primary so- 
lution Zphot ~ 3.10. This is an example of an object with 
significant degeneracy in redshift space: the corresponding 
redshift distribution also shows non-negligible probabilities 
for this object to be either at Zphot ~ 1.8 or at Zphot > 5. 



4 RESULTS 

4.1 The redshift distribution 

Figure |21 shows our derived redshift distribution for Roche 
et al. sample of ERGs. A few objects have been excluded 
either because their light is contaminated by bright neigh- 
bours or because they are quite likely to be stars given their 
very blue (J — Kg) colours. These objects are characterised 
by very low probabilities of being at any redshift as deduced 
by 'hyperz'. We constructed the histogram in fig. EJa) tak- 
ing into account 'hyperz' primary solutions only, for ease of 
comparison with the Hubble diagrams shown in Section 4.2. 
For the ERGs with identification numbers el504 and el605 
in the Roche et al. sample, we adopted the redshifts esti- 
mated by the BPZ code by Bem'tez (2000) instead of the 
'hyperz' primary solution for reasons which are explained in 
Section 4.6. This histogram, as well as all the results pre- 
sented hereafter, include the revised redshifts of the poten- 
tial high-z candidate ERGs e778, elll3 and el272, as they 
were determined after measuring the photometry of these 
objects on the NICMOS and ACS ultra deep images of the 
Hubble UDF. We also show in fig. I^a) the redshift distri- 
bution for the shallower sample of ERGs in the HST-ACS 
field of the GOODS/CDFS obtained by Moustakas et al. 
(2004) (shaded histogram), after division by a factor of 3.08 
in order to normalise their distribution to produce the same 
number of objects in the redshift bin centred at Zphot = 1 as 
found in the present study. 

The redshift distribution for the Roche et al. sample 
spans the range Zphot ~ 0.5 — 4.75. Moustakas et al. sam- 
ple, which is approximately two magnitudes shallower, only 
includes objects up to redshift Zphot = 2.5. The maximum 
of our redshift distribution of Kg < 22 ERGs is located at 
redshift Zphot ~ 1.5. We find a secondary maximum at red- 
shift Zphot ~ 3. However, this secondary peak is not sta- 
tistically significant and disappears when probability densi- 
ties are used to construct the redshift distribution as shown 
in fig. EJb). This latter histogram also includes very small 
corrections for the incompleteness of the sample^ and, thus, 
should be taken as a more realistic representation of the red- 
shift distribution of the ERGs in the GOODS/CDFS deep 
ISAAC field. 



^ Roche et al. ERG sample is considered to be ~ 100% and ^ 80% 
complete to Kg 21.5 and at 21.5 < Kg < 22.0, respectively. 



4.2 The Hubble diagram: Kg vs. Zphot 

Figure |1| shows the Hubble diagram (Kg vs. Zphot) for the 
ERGs in the Roche et al. sample, corresponding to the red- 
shift distribution presented in fig. 1^1 The plots labelled as a) 
and b) show the raw and dust-corrected Kg magnitudes as 
a function of redshift, respectively. Total Kg magnitudes are 
considered in this case. We computed the dust-corrected Kg 
magnitudes using the V-band extinction value Ay applied 
to the best-fit SED in 'hyperz'. The rest-frame wavelength 
sampled by the Kg-band is 



where Aks, is the effective wavelength of the Kg filter, i.e. 

= 2.16 /im. In a dust-screen model, the extinction at 
a wavelength Arf is related to the extinction in the V-band 
(Av) by 

k(Arf) Ay 

^ Ry • (2) 

For the Calzetti et al. (2000) reddening law, Ry = 4.05±0.80 
and k(Arf) is a power law in Arf (cf. 'hyperz' user's manual). 
The dust-corrected Kg magnitude for each source is given 
by the difference (Kg — Aa^j ). It is important to note that 
we assume the validity of the dust-screen model and the 
Calzetti et al. (2000) reddening law for modelling the ex- 
tinction observed in ERGs. We explored the use of other 
reddening laws in 'hyperz' and found that both Milky Way 
and Large Magellanic Cloud-type laws produce very simi- 
lar raw and dust-corrected Hubble diagrams for the ERGs 
in Roche et al. sample. The study of other geometries for 
the distribution of dust is beyond the scope of this paper, 
but one should be aware that they could change the relation 
given by eq. ^ between the extinction Ax^^ and the extinc- 
tion in the V-band Ay. This might have some impact on the 
derived properties of the most highly extincted ERGs. 

In both plots 01 a) and01-b), the circles indicate 'hy- 
perz' primary solutions. The error bars correspond to la 
confidence levels. The overlaid solid curve shows the empir- 
ical K-z relation fitted for massive radiogalaxies (Willott et 
al. 2003), which approximately corresponds to the locus of 
passively evolving present-day 3L* galaxies formed instan- 
taneously at redshift Zf = 10. The dotted line is a nominal 
extrapolation of the same law. In fig. EJa), we have used a 
dashed line to indicate the region of the Hubble diagram 
to which the surveys limited at Kg = 20 have access. Up 
to such a magnitude, only sources with Zphot < 3 are de- 
tected within our survey. For a limiting magnitude Kg = 22, 
the ERGs span the redshift interval Zphot ^ 0.5 - 4.75 and 
display a large dispersion in the Kg vs. Zphot relation. In 
the raw Kg-magnitude Hubble diagram there is an obvious 
lack of objects near the radio-galaxy locus beyond redshift 
Zphot ~ 3. However, comparison with fig. [Ij-b) shows that 
this effect is a consequence of the presence of dust, rather 
than indicating a real absence of intrinsically bright, high- 
redshift ERGs. 

The position of a galaxy in the dust-corrected Kg — Zphot 
diagram depends, of course, on its age as well as its mass. 
The radio-galaxy K-z relation can be interpreted as indicat- 
ing the behaviour of the highest mass galaxies formed at very 
high redshift. However, galaxies may appear to be as bright 
or brighter than high-redshift radio galaxies without being 
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Figure 4. Hubble diagram (total Ka magnitude versus photometric redshifts) for the ERGs in the Roche et al. sample, a) Original Ks 
magnitudes; b) dust-corrected Kg magnitudes. The circles correspond to 'hyperz' primary solutions. The error bars indicate Icr confidence 
levels. The dashed line on a) delimits the region of the Hubble diagram to which Kg = 20 surveys have access. The filled curve corresponds 
to the empirical K-z relation for massive radio galaxies obtained by Willott et al. (2003), which approximately corresponds to the passive 
evolution of a 3L* starburst formed at redshift Zf = 10. The dotted line is a nominal extrapolation of the same law. 



as massive provided their stellar populations are sufficiently 
young and bright. However, as we show in Section 4.5, most 
of the galaxies around the radio-galaxy line do in fact appear 
to at least be massive enough to be the already assembled 
progenitors of the local L > L* population. The comparison 
of fig. EJa) and \^h) also shows that the dispersion in the 
Ks — Zphot relation for ERGs is not an effect produced by 
the presence of different amounts of dust. We see an even 
larger dispersion for the ERGs in the dust-corrected Hub- 
ble diagram. This fact indicates that the ERG population 
comprises objects spanning a wide range in mass. 

To obtain an estimate of the minimum mass of each 
ERG, we used the (k -I- e)-corrected Ks-band absolute mag- 
nitude of each source and computed a lower limit for the lu- 
minosity the galaxy would have at redshift z = 0, assuming 
passive evolution. We explain the details of the calculation of 
the (k-l-e)-corrected Ks-band absolute magnitudes in Section 
4.5. The luminosity L and the absolute magnitude M in a 
given band are related by the relation L/L* = ;^o-°"('^-^*), 
where we considered MJl^j, = —24.2 for h=0.7, from the 2dF 
local luminosity function (Cole et al. 2001). Although the ex- 
act mass-to-light ratio depends on the assumed initial mass 
function (IMF) and the age of each galaxy, we estimate the 
mass of each ERG as (L/L*) x 10^^ M©, where L is the lu- 
minosity the galaxy would have today at redshift z = 0. 
Estimating the masses in this way, we find that the ERGs 
in the Roche et al. sample span two decades in mass, from 
~ 3 X IO^Mq to - 3 X 1O"M0. 



4.3 The red envelope of galaxy evolution 

Figures |3a) and EJb) show the (J — Kg) and (I775 — Ks) 
colours, respectively, as a function of photometric redshift 
for the ERGs in the Roche et al. sample. Only 'hyperz' 
primary solutions are considered in this case. In fig. EJa), 
the filled circles correspond to those objects with J-band 
2"-diameter aperture magnitude above the 2(j-confidence 
limit, in this case J < Jum. = 23.5 (Vega). The up-triangles 
indicate lower limits to the colours of those objects with 
J > 23.5, computed as (23.5 — Ks). In fig.|5Jb), the mixture 
of techniques we applied to measure aperture magnitudes 
on ACS images, i.e. 'SExtractor' plus IRAF 'phot', does not 
allow us to set the 2a confidence limit at a given 1775-band 
magnitude. Thus, we considered as precise colours (filled cir- 
cles) those with 1775-band magnitude error e < 0.5. For ob- 
jects with 1775-band magnitude error e > 0.5 we computed 
lower limits for the colours as (I775 — e — Ks) (up-triangles). 
For comparison, we have added in both plots, |S}a) and|S}b), 
the colours of the three massive galaxies in the HDFS at 
redshift Zphot > 2.4, reported by Saracco et al. (2004) (open 
circles) . 

The extremely red colours observed in ERGs are due to 
two coupled factors: age and dust. As our aim is to study 
which evolutionary model produces the reddest colours ob- 
served in the ERGs, we need to correct for the dust effects 
in order to separate the contribution to the red colours of 
the evolved stellar populations. Figures |H}a) and|5}b) show 
the dust-corrected (J — Ks) and (I775 — Ks) colours, respec- 
tively, as a function of redshift. We computed the dust- 
corrections for the colours subtracting from each magnitude 
the corresponding extinction value given by eq. The 
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Figure 5. (J — Kg) (a) and (I775 — Ks) colours (b) vs. photometric redshifts. The filled circles correspond to exact values and the filled 
up-triangles, to lower limits for the sources with error e > 0.5 in the J or I775 aperture magnitudes. The open circles indicate the colours 
of the three high redshift massive galaxies in the HDFS reported by Saracco et al. (2004). 





Figure 6. Dust-corrected (J — Ks) (a) and (I775 — Ks) colours (b) vs. photometric redshifts. The labels for the filled circles, filled 
up-triangles and open circles are the same as in fig. |F] The different lines show the expected colours for starbursts at different redshifts 
with passive evolution thereafter and metallicity Z = Zq: Zf = 5 (solid), Zf = 4 (dashed-dotted) and Zf = 3 (dotted) . The dashed line 
corresponds to a similar starburst at redshift Zf = 5, but with a higher metallicity Z = 2.5 Z©. 



symbol labels in both plots of fig. |^ are equivalent to those 
in fig. We used the public code 'GALAXEV (Bruzual 
and Chariot 2003) to compute the expected (J — Ks) and 
(I775 — Ks) colours of different kinds of dust- free passively 
evolving galaxies, in order to obtain a suitable model for 
the red envelope delimited by the ERG largest colours. We 



found that an instantaneous starburst at redshift Zf = 5 
with passive evolution thereafter and solar metallicity (solid 
lines in models quite well the dust-corrected red en- 

velope. A similar starburst formed at redshift Zf = 5, but 
with a higher metallicity Z — 2.5 Zq (dashed line), has been 
added for comparison. The plots in fig. |H] show that so- 



© 0000 RAS, MNRAS 000, 000-000 



8 K.I. Caputi, J.S. Dunlop, R.J. McLure and N.D. Roche 



lar metallicity Z = Z0 is sufficient to explain the reddest 
colours observed. Even the few apparent outliers at redshift 
Zphot ~ 0.5 — 2.0 can be explained by such a model once 
the error bars are taken into account, which are ~ 0.2 mag 
for these objects. On the other hand, the modelling of the 
expected colours for starbursts at different redshifts allows 
us to investigate the epoch since when the ERGs lying at 
redshifts Zphot > 2 have been following passive evolution. 
The dotted and dashed-dotted lines in fig. |S] correspond to 
the modelled colours for starbursts at redshifts Zf = 3 and 
zf = 4, respectively (both with passive evolution thereafter 
and solar metallicity Z — Zq). We see that the reddest ERGs 
lying at redshifts Zphot ~ 2 (3) have been passively evolving 
since redshift Zphot ~ 3 (4). This does not necessarily mean 
that these objects have been formed in these epochs, but 
corresponds to the formation redshifts of the youngest stel- 
lar population present in them. 

The comparison of fig.|S]and fig.|S]provides a simple and 
useful way to understand the extent to which the presence 
of dust is responsible for the original extremely red colours 
of ERGs. This fact is summarised in fig. |71 where we show 
the median of the V-band extinction values Ay necessary 
to deredden the ERG best-fit SEDs at different redshifts. 
We see that considerable median extinction values (Av > 
0.5) are necessary to explain the red colours observed in 
ERGs up to redshift Zphot ~ 4. For objects below redshift 
Zphot = 1, the modelled Ay values are very large, illustrating 
that only extremely dusty starbursts can classify as ERGs at 
relatively low redshifts. The median of the extinction Ay is 
lower at redshifts Zphot 1 — 2, where perhaps the essentially 
non-dusty evolved stellar populations are expected to make 
the most of their contribution. At redshifts Zphot ~ 2 — 4, 
where according to fig. Inmost of the objects seem to have 
experienced recent starbursts, the median of the extinction 
Ay increases again. For objects beyond redshift 

Zphot — ^, 

nearly no extinction is required indicating that, at this stage, 
the extremely red colours are mainly due to the shift of the 
Lyman break to optical wavelengths rather than dust. 

4.4 The evolving luminosity function of ERGs 

Figures |5] and |5] show the luminosity function (LF) of the 
ERGs in Roche et al. sample in different redshift intervals. 
We assigned to each object the Kg-band k-corrected absolute 
magnitude (Mks) obtained as part of the 'hyperz' output. 
We considered 'hyperz' primary solutions for binning the 
data in redshift space. We computed the comoving densi- 
ties per absolute magnitude bin by dividing the number of 
sources in each bin by the corresponding comoving volume. 
To take into account the limits of the survey (Kg — 22), 
we corrected the contribution of each source by a weighting 
factor Vmaxbin/Vmaxobs, whcrc Vmaxbin is the volumc deter- 
mined by the maximum redshift of the bin and Vmaxobs is 
the volume corresponding to the maximum redshift at which 
the source would still be included in the survey (provided 
it is lower than the maximum redshift of the bin). We also 
applied a correction factor to account for a slight incom- 
pleteness of the sample at 21.5 < Kg < 22.0. 

Fig. m shows the ERG LF at redshifts (zphot) = 1.0, 
(zphot) = 1-5 and (zphot) ~ 2.0, indicated by circles, squares 
and up-triangles, respectively. The error bars correspond to 
the maximum of the Poissonian errors and the errors due 




Figure 7. Median of the V-band extinction values applied to the 
best-fit SEDs vs. photometric redshifts. Only 'hyperz' primary so- 
lutions are taken into account. The horizontal and vertical error 
bars indicate the binning in redshift space and the largest bin- 
ning used in 'hyperz' for the iteration through the Ay values, 
respectively. 



to cosmic variance, which we take on average as 40% in 
the number counts at these redshifts (cf. Somerville et al. 
2003). The inclusion of the cosmic variance is fundamen- 
tal for highly-clustered populations, as ERGs at redshifts 
z ~ 1 — 2 are known to be. The global Kg-band LF at 
redshifts 1 < Zphot < 2 has been added for comparison. 
The solid line represents the average Kg-band LF of the 
HDFS and HDFN, obtained by integrating in each abso- 
lute magnitude bin the best Schechter function fitted by 
Bolzonella, Pello & Maccagni (2002) for datasets in these 
fields up to a limiting magnitude of Kg ~ 23. The dotted 
lines indicate the largest error bars in the normalisation 
parameter 0*, corresponding to the fitting of the HDFS 
data. As expected, our ERG LF reproduces the shape of 
the bright end of the global Kg-band LF. The differences 
at fainter magnitudes are at least in part due to the na- 
ture of the ERG colour selection. We find only 3 very lu- 
minous galaxies with —27 < Mk^ < —26 with estimated 
redshifts Zphot e [0.75;2.25], 2 of which have redshifts in the 
interval [1.75;2.25]. These 2 sources have (I775 - Kg) > 5.9 
and (J — Kg) > 2.3. For the absolute magnitude range 
—26 < < —24, we do not observe any evolution in 

the ERG LF from redshifts ( 

Zphot 

) = 1.0 to (z 

phot ) = 2.0. 

Any fluctuation in the mean values of the LF in this abso- 
lute magnitude bin can be accounted for within the cosmic 
variance error bars. For absolute magnitudes —24 < Mk^ 
we do not observe any evolution in the LF from redshifts 
(zphot) ~ 1-0 to (zphot) ~ 1.5. The limits of the survey do 
not allow us to explore the LF at redshift (zphot) ~ 2.0 for 
such faint objects. 

In figure |5| we show again the ERG LF at redshifts 
(zphot) = 1.0, 1.5, 2.0 and we add the ERG LF at (zphot) = 
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Figure 8. Luminosity function of ERGs at different redshifts: 
(zphot) = 1-0 (circles), (zphot) = 1-5 (squares) and (zphot) = 2.0 
(up-triangles). The average global Kg-band LF in the HDFN and 
HDFS has been overlaid for comparison (solid line). The latter 
has been computed integrating in each absolute magnitude bin 
the Schechter functions fitted by Bolzonella et al. (2002) in these 
fields. The error bars have been estimated using the maximum 
uncertainty in the normalisation parameter </)*, obtained for the 
HDFS (dotted lines). 



2.5, for comparison. We computed the latter taking into ac- 
count all the objects in the redshift interval [2.0,3.0]. At 
these redshifts, we have adopted a slightly lower value for the 
typical cosmic variance, only 30% in the number of counts 
on average (Somerville et al. 2003). This plot shows that 
the bright end of the ERG LF does not decrease significantly 
from redshifts (zphot) ~ 2.0 to (zphot) = 2.5. This result con- 
firms the existence of a population of extremely red bright 
galaxies at high redshifts. We find 7 objects with estimated 
redshifts in the interval [2.0,3.0] and absolute magnitudes 
-27 < Mk, < -26, 6 of which have (I775 - K^) > 4.8, plac- 
ing them among the very reddest objects at these redshifts. 
Moreover, all of these objects have (J — Ks) > 2.2. In the 
next section we investigate the change of the absolute mag- 
nitudes of ERGs after evolution to redshift z = and explore 
the comoving densities of potential progenitors of the local 
L > L* population. 



4.5 Comoving densities of massive galaxies 

In this section we study the comoving densities of ERGs 
that will still contribute to the bright end of the Ks-band 
LF once they evolve down to redshift z = 0. If we assume 
passive evolution for ERGs from the epoch of observation, 
the maximum dimming of the absolute magnitude Mk^ can 
be estimated and, thus, a lower limit to the luminosity each 
object would have at redshift z=0 obtained. To do this, we 
need to estimate an evolutionary correction (e-correction) 
factor at redshift z = which, when added to the 'hyperz' 



Figure 9. Luminosity function of ERGs at different redshifts. 
The references are the same as in 

fig- m The LF at redshifts 
2 < Zpijot < 3 has been added (star-like symbols) . This plot shows 
that the bright end of the ERG LF does not obviously decrease 
from redshifts (zphot) = 2.0 to (zpijot) = 2.5. 



k-corrected absolute magnitude, yields the absolute mag- 
nitude the galaxy would have today. The e-correction for 
each galaxy depends on its spectral type and age. We used 
the public code 'GALAXEV (Bruzual and Chariot 2003) to 
compute the absolute magnitudes at redshift z = of the 
spectral types corresponding to our ERGs at different ages. 
The correction applied to each object is the difference in 
the absolute magnitude the object would have today and at 
the age it had when its light was emitted. The objects that 
are considered to be progenitors of the local L > L* popu- 
lation are those with (k -|- e)-corrected absolute magnitude 
Mk, < Mk^, given the relationship L/L* = ;^o-° *(^-^*'. 
In a passive evolution scenario, these ERGs must have con- 
tained a minimum stellar mass of ~ 10^^ Mq when their 
light was emitted. 

Figure lTUl shows the comoving densities at different red- 
shifts of the 27 ERGs in Roche et al. sample which are ex- 
pected to have L > L* at redshift z = under passive evo- 
lution (filled circles). We used M^^^ = -24.2 for h=0.7 (Cole 
et al. 2001). Once more, we computed the comoving densi- 
ties dividing the number of sources in each redshift bin by 
the corresponding comoving volume. We divided the sam- 
ple in bins of width Azphot = 0.5 up to redshift Zphot ~ 2 
and Azphot ~ 1 for redshifts 2 < Zphot < 4. To compute 
these comoving densities, we did not use 'hyperz' primary 
solutions but a redshift probability density distribution that 
we constructed normalising for each source the percentage 
probabilities given in the 'hyperz' output file Jog. phot. In 
this way, each source may contribute in a fraction to different 
redshift bins. Only sources with 'hyperz' primary solutions 
Zphot < 4 have been taken into account. We applied a weight- 
ing factor of the kind Vmaxbin/Vmaxobs, as it was explained 
in Section 4.4, to the normalised probability density distri- 
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Figure 10. Lower limits on the comoving densities of ERGs 
progenitors of local L > L* galaxies as a function of redshift 
(filled circles). The horizontal error bars indicate the binning in 
the redshift space. The vertical error bars correspond to Pois- 
sonian errors in the number of objects. The lower limit on the 
comoving density of progenitors of L > L* galaxies at redshift 
(zphot) = 2.7 obtained by Saracco et al. (2004) (open circle) and 
the comoving density of local L > L* galaxies (star-like sym- 
bol) have been added for comparison. The dashed-dotted and 
dashed lines show the comoving densities of dark matter halos 
with masses M > 2.5 X 10^2 and M > 1.0 X lO" Mq, respec- 
tively, as obtained from ACDM models of structure formation. 



bution of each source. However, in each case, the weighting 
factor has been estimated using only the k-corrected ab- 
solute magnitude Mk^ given by 'hyperz' in the output. A 
more rigurous procedure should take into account a prob- 
ability density distribution for the absolute magnitudes in 
correspondence with the probability density distribution for 
redshifts. We also corrected for the incompleteness of the 
sample. The error bars for the comoving densities corre- 
spond to Poissonian errors in the number of sources, which 
are dominant in this case. We added for comparison the co- 
moving density of the local L > L* population (star-like 
symbol), which we computed as the average of the values 
obtained integrating the local Ks-band LFs fitted by Cole et 
al. (2001) and Kochanek et al. (2001). 

Our results show that there is no significant change in 
the comoving densities of ERG progenitors of L > L* galax- 
ies between redshifts 1.5 < Zphot < 2.0 and 2 < Zphot < 3, 
which is consistent with the results presented in Section 4.3 
for the bright end of the LF. Actually, all the galaxies with 
—27 < < —26 in Section 4.3 are included in the progen- 
itors of L > L* (z = 0) subsample. The values we obtain for 
the comoving densities are: Pc = (8.5±3.5) X 10"^ Mpc-3 for 
(zphot) = 1.25, pc = (6.1 ± 2.8) X 10-5 Mpc-^ for (zphot) = 
1.75, pc = (6.1 ± 1.9) X 10-5 Mpc"^ for (zphot) = 2.5 and 
Pc = (2.1 ± 1.1) X 10-5 Mpc"^ for (zphot) = 3.5. It is im- 
portant to note that the relatively small comoving densities 
at (zphot) ~ 1-25 and (zphot) ~ 1-75 (more than 10 times 



smaller than the corresponding value we adopted for local 
L > L* galaxies) could in part be due to the ERG colour 
cutoff. If this is the case, we should conclude that the ERGs 
at redshifts z ~ 1 — 2 cannot account for the whole popula- 
tion of progenitors of the local massive galaxies. Although 
the ERGs reproduce quite well the shape of the bright end of 
the global Ks-band LF, the values are lower, indicating that 
larger comoving densities of massive objects are expected to 
be obtained in a Kg-selected sample without a colour selec- 
tion. Thus, the conclusion of a constant comoving density 
of massive objects from redshifts 1.5 < Zphot < 2.0 and 
2 < Zphot < 3 is in principle only applicable to ERGs, and 
may of course not hold for a general Kg-selected sample. The 
mean value we obtain for the lower limit to the comoving 
density of L > L* progenitors at (zphot) ~ 2.5 is only one 
half of the mean value obtained by Saracco et al. (2004) at a 
similar redshift (open circle in fig. llOH . but there is no signif- 
icant discrepancy between the results given the large error 
bars. We conclude that Saracco et al. large mean comoving 
density might be due to a simple cosmic variance effect and 
much larger samples of Ks-selected galaxies are necessary to 
establish a robust estimate of the fraction of the stellar mass 
already assembled in ERGs at these redshifts. ^ 

In fig. llOl we show the comoving densities of dark matter 
halos with masses M > 1.0 X 10" Mq (dashed lines) at dif- 
ferent redshifts, as they are predicted by ACDM models of 
structure formation (Kauffmann & Chariot 1998, Somerville 
& Primack 1999). We also show the comoving densities of 
dark matter halos with masses M > 2.5 X 10^^ Mq (dashed- 
dotted lines), a mass threshold deliberately selected to co- 
incide with the comoving densities of L > L* galaxies at 
z = 0. At any redshift, our comoving densities of progenitors 
of local massive galaxies are lower than the corresponding 
comoving densities of dark matter halos massive enough to 
host these galaxies. For instance, the comoving densities of 
M > 2.5 X 10^^ Mq halos decreases less than a factor ~ 2 
from redshifts z = 0toz=1.5, while the comoving densi- 
ties of ERGs progenitors of local massive galaxies are more 
than 5 times below this value. This leads to two extreme 
interpretations: i) even taking into account the possible pro- 
duction of mergers (i.e. the increasing comoving densities of 
dark matter halos more massive than a given mass thresh- 
old with time) , ERGs do not seem to be sufficient to account 
for all the progenitors of local L > L* galaxies. This would 
reinforce the idea that the ERG colour cutoff must be at 
least partially responsible for this deficiency, which appears 
to be more important at lower (zphot ~ 1 — 2) than higher 
(zphot ~ 3 — 4) redshifts. ii) ERGs do account for the pro- 
genitors of local massive galaxies but the passive evolution 
scenario is not completely valid. Galaxies which are not suf- 
ficiently massive at a given redshift might build up more 
stellar mass at later epochs and, thus, should also be con- 
sidered as progenitors of the local L > L* galaxy population. 
The massive ERGs at z > 2 could instead be progenitors of 

Saracco et al. claim to reproduce about 40% of the comoving 
density of the local massive galaxies, after comparison with the 
local comoving density of L > L* early-type galaxies obtained 
from the LF fitted by Marzke et al. (1998) using B-band data. 
Instead, we prefer to compare the comoving densities with the 
local values obtained integrating more recent LFs directly fitted 
on K-band data (Cole et al. 2001, Kochanek et al. 2001). 
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much more massive (~ 4L*) local galaxies, in which case our 
results show that all of them are already in place at redshift 
z ~ 3.5. 

Interestingly, the progenitors of L > L* galaxies are 
among the reddest ERGs. At different redshifts, these mas- 
sive galaxies in the Roche et al. sample have ((I775 — Ks)) ?s 
4.8 — 5.1. The mean (J — Kg) colours increase with redshift, 
as it is expected when the 4000A break enters the region be- 
tween the J and Kg filters. In all cases, (J — Kg) > 2. Some 
of these massive galaxies are among the reddest objects in 
(J — Kg) at redshifts Zphot > 2, due to the combined effects 
of the break and considerable amounts of dust (wc find V- 
band extinction values of up to Ay = 1.8 for these mgissive 
objects). The characteristics of these sources are consistent 
with the 'hyper extremely red object' population discovered 
by Totani et al. (2001). However, not all the (J — Kg) red- 
dest objects evolve to L > L* sources in a passive evolution 
scenario. In some cases, the extremely red (J — Ks) colours 
are almost exclusively due to very large dust extinctions. 
Recently, some extremely red (J — Ks) sources have been 
found to be counterparts of sub-millimetre sources, which 
are known to be mainly located at redshifts 2 < z < 4 
(Dunlop 2001, Aretxaga et al. 2003, Chapman et al. 2003). 
However it is not yet clear whether the ERGs associated 
with sub-millimetre sources are primarily high or low mass 
objects (Prayer et al. 2004). 

4.6 Sources beyond redshift z = 4 - the ERG - 
Lyman-break connection 

In this section we focus on the properties of the Lyman 
break galaxies in the GOODS/CDFS included in the Roche 

et al. ERG sample. These objects could be candidates for 
massive galaxies at very high redshift. liere we summarise 
the detailed properties of the individual objects and discuss 
the probability that they are genuinely located at such high 
redshifts. 

Several studies have already been published on the 
selection of very high redshift candidate sources in the 
GOODS fields (Stanway, Bunker & McMahon 2003; 
Bouwens et al. 2003; Bremer et al. 2004; Dickinson et al. 
2004) . In most of these studies a colour cutoff selection tech- 
nique has been applied (usually (I775 — z) > 1.3 — 1.5, AB) 
to select candidates at redshifts z > 5. We find 9 of these 
candidate sources in the Roche et al. ERG sample. How- 
ever, we cannot confirm any of these sources as being at 
very high redshift. We argue that the fraction of low- mass 
star and lower-redshift ERG contaminants is usually under- 
estimated when a single-colour selection criterion is applied. 
Multicolour photometry appears to be a much more power- 
ful way to obtain reliable candidate lists of very high redshift 
objects, and allows the investigation of degerate solutions in 
redshift space. 

Using 5-epoch GOODS HST-ACS (Vl.O) and VLT- 
ISAAC data, we obtained an initial 'hyperz' output list con- 
taining 12 sources with a primary solution at Zphot > 4, 
which adopted as candidates for detailed inspection. We 
present the properties of the complete list of 12 candidate 
sources in Table 1. We used the public code BPZ (Bem'tez 
2000) to obtain estimated redshifts for these objects in an in- 
dependent way. We also made an individual study of each of 
our 12 candidate sources, inspecting their ACS images and 



analysing their magnitudes and colours in different bands. 

After a first analysis, we rejected the following 7/12 z > 4 
candidates for the reasons explained below: 

• ell4-e566-e967-e2386 : these sources appear to be 
at redshifts Zphot > 5.3 from 'hyperz' and all of them have 
(I775 — z) ~ 1.3 — 1.5, AB. Only e566 is not confirmed by the 
BPZ code as a very high redshift source (it is estimated to 
be at redshift Zphot = 1.28). However, the probabilities as- 
sociated with the 'hyperz' primary solutions, obtained from 
the minimisation, are very low (P < 19%) in all cases. 
e967 corresponds to source number 4 in Stanway et al. list 
and source number 8 in the list of Bremer et al. el 14 and 
e2386 are sources number 13 and 9 in the Bremer et al. 
list, respectively. These sources are unresolved on ACS im- 
ages and have quite blue near-IR colours, (J — Ks) < 1.7 
(Vega) in all cases. e967 and e2386 are significant detec- 
tions in the V-band (V = 27.9 ± 0.4 and V = 27.4 ± 0.3, 
AB, respectively) and e566, in the B-band (B = 23.7 ± 0.1, 
AB). Thus, it is more likely that they are cool staxs than 
high redshift galaxies. Indeed, e967 has been spectroscop- 
ically confirmed as such (cf. Dickinson et al. 2004). These 
cool stars are known to be among the main contaminants of 
very high redshift candidates. 

• el504: this source has 'hyperz' primary solution 
Zphot — 4.63, but actually has a probability P s at any 
redshift between and 10. The BPZ code gives an estimated 
redshift of Zphot = 1-44 for this source. Besides, although it 
is not detected in the V-band it is a significant detection 
in the B-band (B = 27.2±0.3, AB). This source corresponds 
to object number 9 in the Stanway et al. list. 

• el605: this source has 'hyperz' primary solution 
Zphot = 6.71 (with an associated probability P ^ 90%), but 
the BPZ code gives Zphot = 2.84. This source is not detected 
on the z-band images, but it could have some flux in the I776- 
band (I775 = 27.9±1.3, AB). Its very red near-IR colour, i.e. 
(J — Ks) > 2.5 (Vega), suggests that the estimated redshift 
from the BPZ code is more reliable in this case. 

• e2006: this source has two 'hyperz' solutions with sim- 
ilar probabilities. The primary solution is Zphot — 6.21 (P « 
73%) and the secondary one is Zphot = 1.15 (P ^ 60%). The 
BPZ estimation is Zphot = 6.37. We measured (I775 — z) = 
2.1 ± 0.8 for this object and we do not find it is significantly 
detected in either the B or V bands. However, the higher red- 
shift solution would imply it is an extremely bright galaxy 
with a k-corrected absolute magnitude Mk^ = —28.5 and 
estimated mass ~ 6.7 X IO^^Mq, so we consider that in this 
case the 'hyperz' secondary solution is more likely to be the 
right answer. 

3/12 high-z candidate ERGs lie in the Hubble UDF, which 

images have been recently made publicly available. After 
revision of the photometry on the HST-ACS and HST- 
NICMOS ultra-deep optical and near-IR images, we reject 
also these 3/12 high-z candidates: 

• e778: using 5-epoch GOODS HST-ACS (Vl.O) and 
VLT-ISAAC data, the 'hyperz' primary solution for this 

source is Zphot = 4.74 (with an associated probability P ~ 
98%). The BPZ redshift estimation is Zphot = 5.42. However, 

^ even when its flux is manually measured with the IRAF task 
'phot' 
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Table 1. List of ERGs in the Roche et al. sample with 'hyperz' primary solutions Zpjiot > 4. The first column is the identification 
number. The second and third columns indicate the coordinates of the source obtained from the Ks-band images. The fourth and fifth 
columns show the redshifts estimated by 'hyperz' and the BPZ code, respectively. Column six lists comments on the individual objects. 



ERG id RA(J2000) DEC(J2000) hyperz BPZ accepted? COMMENTS 



114 


3 


32 


22 


46 


-27:50 


47 


16 


5 


gQ+0.03 
°'^-0.14 


5 


75 


N 


Unresolved. (J - Ks)vcga = 0.8 ± 0.3 


566 


3 


32 


24 


78 


-27:49 


12 


91 


6 


gg + 0.11 
"='-0.12 


1 


28 


N 


Unresolved. (J - Ks)voga = 1.1 ± 0.3. Bab = 23.7 ± 0.1 


967 


3 


32 


18 


17 


-27:47 


46 


48 


5 


on+0.03 

•^^'-o.os 


5 


70 


N 


Spectroscopically confirmed star. Vab = 27.9 ± 0.4 


1504 


3 


32 


18 


17 


-27:46 


16 


33 


4 


cQ+0.10 
"'-'-0.06 


1 


44 


N 


Bab = 27.2 ±0.3 


1605 


3 


32 


19 


67 


-27:46 


02 


04 


6 


7-1+0.75 
' ^-0.32 


2 


84 


N 


(J — Ka)vega > 2.5. Not detected in z band. 


2386 


3 


32 


19 


23 


-27:45 


45 


39 


5 


OO+0.05 
•-''^-0.04 


5 


56 


N 


Unresolved. (J - Ks)vcga = 1.2 ± 0.2 


2006 


3 


32 


28 


81 


-27:44 


30 


54 


6 


91+0.17 
^^-0.33 


6 


37 


N 


Significant secondary solution Zpi^ot = 1.15 


778 


3 


32 


41 


76 


-27:48 


24 


92 


4 


74+0.20 
'^-0.25 


5 


42 


N 


Revised UDF photometry implies Zpi^ot = 1.56 


1113 


3 


32 


34 


65 


-27:47 


20 


89 


5 


07+0.16 
-0.08 


1 


05 


N 


Revised UDF photometry implies Zp^ot = 1.43 


1272 


3 


32 


41 


69 


-27:46 


55 


40 


5 


nQ+0.21 
"'-'-0.13 


5 


48 


N 


Revised UDF photometry implies Zpjjot = 1.91 


491 


3 


32 


25 


95 


-27:49 


30 


38 


4 


c-o+0.24 
^°-0.09 


5 


22 


Y 


No significant 'hyperz' secondary solutions. 


1703 


3 


32 


37 


74 


-27:45 


05 


41 


4 


7^+0.04 
'^-0.05 


5 


38 


Y 


No significant 'hyperz' secondary solutions. 




Figure 11. Redshift probability density distributions for the 
accepted z > 4 candidate sources in the Roche et al. ERG sample. 



This object has (I775 — z) = 0.6 ± 0.3 and it is not detected 
in either the B or V bands. 

• el703: this source has a 'hyperz' primary solution 
Zphot ~ 4.74 ( P ~ 90%) and no significant secondary so- 
lution. The BPZ redshift is Zphot ~ 5.38. The estimated 
Kg-band k-corrected absolute magnitude for this source is 
= —27.3 and the estimated mass is 2.7 X 1O"M0 . It 
has (I775 — z) = 0.97 ± 0.14 and it is not detected in either 
the B or V bands. 

Finally, we list other ERGs in the Roche et al. sample 
which are among the very-high redshift candidates selected 
by other authors, and for which we obtain lower redshift 
estimates: 



after revision of the photometry on the Hubble UDF images, 
this source appears to be a significant detection in the B- 
band. The revised 'hyperz' primary solution is Zphot ~ 1.56. 

• elll3: this source initially had a 'hyperz' primary so- 
lution Zphot = 5.07 (P ~ 56%) and secondary solution 
Zphot = 1.10 (P « 26%). The BPZ code favoured a low 
redshift solution Zphot = 1.05. We did not detect this source 
in either the B or V bands on the GOODS VI. images. 
However, it is a clear detection in both the B and V bands 
on the Hubble UDF images. The revised 'hyperz' primary 
solution is Zphot = 1.43. This source corresponds to object 
number 2 in Stanway et al. list. 

• el272: this source initially had a 'hyperz' primary so- 
lution Zphot = 5.03 ( P « 82%) and no significant secondary 
solution. The BPZ estimation was Zphot = 5.48. It was not 
a significant detection in either the B or V bands on the 
GOODS VI. images, but it is on the Hubble UDF maps. 
The revised 'hyperz' primary solution is Zphot = 1.91. 

We accept the remaining 2/12 as z > 4 candidates, whose 
redshift probability density distributions we show in fig. 1111 

• e491: 'hyperz' primary solution is Zphot ~ 4.58 (with 
an associated probability P ~ 98%). The BPZ code esti- 
mates a redshift of Zphot = 5.22. There are no significant 
secondary solutions for the redshift of this source in 'hyperz'. 
The Ks-band k-corrected absolute magnitude for this galaxy 
is Mk, = —25.6 and we estimate a mass of 5.8 x lO^'^Mm. 



ERG id 


Author 


Author's id 


'hyperz' 


82 
309 
225 
1423 


Bremer et al. 
Bremer et al. 
Dickinson et al. 
Bremer et al. 


7 
6 

SiDOlS 
2 


1 40+0-30 

J^-^'J-0.04 
9 09+0.11 
^"^^-0.05 

1 74+0.08 

'^-0.10 

1 60+0-28 

l.DU_o Qg 



From the sources listed above, only e309 has a very-high- 
redshift secondary solution, although with negligible signif- 
icance, at Zphot ~ 6.74. 

To summarise, we conclude that, within the limits of 
our survey, only a few ERGs are likely to be at very high 
redshifts. Deeper near-IR selected samples of galaxies are 
necessary to confirm or refute the rarity of massive objects 
assembled at very early epochs. 



4.7 Review of the ERG cluster candidate at 
redshift z ~ 1.5 

Roche et al. (2003) found 10 ERGs within a 20" radius 
of the Chandra source XID:58 (RA=3:32:11.85, DEC=- 
27:46:29.14, J2000), which is itself a counterpart of the ERG 
identified as el435 in Roche et al. list. This overdensity sug- 
gested the presence of an ERG cluster at redshift z ~ 1.5, 
based on the photometric redshift of the Chandra source 
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XID:58, estimated as Zphot = 1-44 by Mainieri (2003). Based 
on the source colours, Roche et al. argued that 7/9 of the 
ERGs around XID:58 (el435) could be part of that clus- 
ter. In the following table, we present our redshift estimates 
(both primary and secondary) for el435 and the 7 surround- 
ing sources that could be associated with it. 



ERG id Zphot (primary) Zphot (sec.) 



1333 


3 10+°" 


1.78 


1341 


1 (-9+0. 06 


1.01 


1390 


1 2C)+"-"i 

^^•^='-0.01 


1.67 


1404 


, (-9+0.02 


1.78 


1423 




1.34 


1435 (XID:58) 


1 or+O.OS 
1.00_Q 08 


1.30 


1464 


o iq+o.07 
o-J^y_o.io 


1.48 


1481 


1 sn+012 

1.0U_Q ^2 


1.01 



Our estimated redshifts seem to confirm the presence of 
a cluster at redshift Zphot ~ 1-5. However, it is not clear 
whether all the sources listed above belong to it. el333 
and cl464 seem to be higher redshift sources, although the 
'hyperz' secondary solution might locate el464 at redshift 
Zphot = 1-48. For el435 (counterpart of XID:58) we esti- 
mate a redshift Zphot = 1-85, a value somewhat higher than 
Mainieri's estimation. Our secondary solution, on the con- 
trary, is closer to their value. The sources which seem to be 
part of a cluster at redshift Zphot ^ 1.5 are el341, el404, 
el423 and el481. 

Roche et al. also suggested that the 2/9 remaining 
ERGs around XID:58 could be associated with another X- 
ray source (Chandra X1D:149), which is spectroscopically 
confirmed to be at redshift Zphot = 1.033 (Szokoly et al. 
2003). Our estimated redshifts for these 2 remaining ERGs 
are: 



ERG id 


Zphot (primary) 


Zphot (sec.) 


1311 
1467 




0.61 

.").21 



The 'hyperz' primary solutions seem to confirm that these 
2 ERGs are actually associated with the X-ray source 
XID:149. 



5 SUMMARY AND CONCLUSIONS 

In this paper we have presented estimated redshifts for the 
Extremely Red Galaxies selected by Roche et al. (2003) in 
the 50.4 arcmin^ of the GOODS/CDFS deep ISAAC field 
to a limiting magnitude Kg = 22. This is the deepest sig- 
nificant sample of ERGs selected to date and constitutes a 
complement to other shallower but wider surveys of ERGs 
in the same field (Moustakas et al. 2004). 

We have used multicolour photometry in seven pass- 
bands (B,V, I775, z, J, H and Kg) to compute photometric 
redshifts using the public code 'hyperz'. The B,V, I775 and 
z magnitudes have been measured on the stacked 5-cpoch 
GOODS ACS images (VI. data release). We have obtained 



the redshift distribution and the corresponding Hubble di- 
agram (Ks — Zphot), which show the existence of ERGs up 
to redshifts Zphot ~ 4.75 at this depth. The ERG popula- 
tion is characterised by a large dispersion in the Ks — Zphot 
relation. We find that this dispersion is even more impor- 
tant in the dust-corrected Hubble diagram, indicating that 
the ERG populaton is composed of objects spanning a wide 
range in mgiss. We estimate that the ERGs in Roche et al. 
sample span two decades in mass, from ~ 3 x IO^Mq to 
~ 3 X 10" Mq. 

We have studied the red envelope of galaxy evolution, 
determining the galaxy template which best describes the 
ERG dust-corrected reddest colours observed in the Roche 
et al. sample as a function of redshift. We find that a star- 
burst formed at redshift Zf = 5 with passive evolution there- 
after and solar mctalhcity provides a very good description 
of the red envelope of ERG evolution for objects selected 
with Ks < 22. Our work indicates that the simplified tra- 
ditional picture for ERGs as either old elliptical galaxies 
or young dusty starbursts can no longer be sustained. The 
colours of many of the reddest ERGs can only be explained 
by the combined effects of evolved stellar populations and 
dust. 

We have also explored the evolution of the ERG lumi- 
nosity function from redshifts (zphot) = 1 to (zphot) = 2.5. 
The ERG LF at redshifts 1 < Zphot < 2 reproduces the 
shape of the bright end of the global Ks-band LF, as per- 
haps expected. We find no evolution in the bright end of the 
LF of ERGs from redshifts (zphot) = 1-0 to (zphot) = 2.5 
(-26 < < -25) and from redshift (zphot) = 2.0 to 

(zphot) = 2.5 (-27 < Mks < -26). This effect cannot be 
deduced from inspection of the Hubble diagram for ERGs, 
constructed using the raw Kg magnitudes. Only after the 
appfication of dust corrections are some ERGs revealed as 
bright and massive galaxies. 

One of the main goals in the selection of ERG samples 
is to trace the progenitors of very luminous local galaxies 
(L > L*). Only a small fraction (~ 15%) of the ERGs in the 
Roche et al. sample can evolve into such luminous galax- 
ies under passive evolution. At redshifts Zphot 1 — 2, we 
find that the comoving density of ERGs progenitors of mas- 
sive galaxies is less than ~ 1/10 of the local value. Between 
redshifts z = and z = 1.5, ACDM models only predict 
a decrease in the comoving densities of dark matter halos 
massive enough as to host local L > L* galaxies by a factor 
~ 1.5-2. This suggests that either ERGs cannot account for 
the whole population of local massive galaxies or the passive 
evolution scenario is not completely valid. This deficiency 
appears as being more dramatic at lower (zphot ~ 1 — 2) 
than higher (zphot ~ 3 — 4) redshifts, suggesting that ERGs 
could be rather better tracing the progenitors of local mas- 
sive galaxies present at early epochs. 

Finally, we have investigated the existence of Lyman 
break gala^xies massive enough to be included in an ERG 
sample. Out of an initial candidate list of 12 sources, we 
only confirm 2 ERGs as having high probabilities of being 
located at redshifts Zphot > 4. Deep observations in wider 
areas arc crucial in order to constrain lower limits on the 
fraction of massive galaxies with the bulk of their stellar 
mass already assembled at very high redshifts. 

Several open questions still remain for the ERG pop- 
ulation in the GOODS/CDFS field. An important issue is 
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to understand how ERG morphology changes with redshift, 
a problem we plan to study in a future paper. Also, very 
deep spectroscopic observations scheduled with the Gemini 
Multiobject Spectrograph (GMOS) will provide more accu- 
rate redshifts for a subset of the ERGs in the Roche et al. 
sample, allowing a detailed study of their properties. 
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